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Abstract. Here we are concerned about uniform stability of damped nonlinear trans- 
verse vibrations of an elastic string fixed at its two ends. The vibrations governed by 
nonlinear integro-differential equation of Kirchoff type, is shown to possess energy uni- 
formly bounded by exponentially decaying function of time. The result is achieved by 
considering an energy-like Lyapunov functional for the system. 
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1. Introduction 

The textbook treatment of transverse vibrations of a metallic wire (string), governed by 
the linear wave equation, does not portray existence of whirling out of plane motion. The 
phenomenon was first observed by Hunton as reported by Harrison [9|. The reference 
to other subsequent experiments can be found in the book by Nayfeh and Mook 1141 . 
It is observed that the whirling motion, sometimes also referred to as ballooning motion 
occurs when the amplitude and frequency of a plane excitation and the phase difference 
between the response and the excitation exceeds certain critical values. The explanation 
of the phenomenon lies in a nonlinear treatment of the problem. Long ago Kirchoff 1 10 1, 
taking longitudinal elasticity also into consideration, gave the following nonlinear integro- 
differential equation for transverse vibrations confined to a plane: 



where u(x,t) is the transverse deflection, I the length of the string, 8 > the coefficient of 
damping, a 2 = To/pA and b = E/2pl. The constants Tq,p,A,E are respectively the initial 
axial tension, the mass density, the cross-section area and Young's modulus of the string. 
For treating nonlinear out of plane motion with components (v, w) in the perpendicular y 
and z directions, Anand [ 1 1 derives a generalization of equation (la): 




< x < I, t > 



(la) 




(lb) 
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where u has components (v,w). A systematic rederivation of (lb) can also be found in 
Nayfeh and Mook 1 14 1. The book [14| also refers to other authors who derived different 
nonlinear differential equations under various assumptions. 
For a string fixed at both ends, we have 

u(0,?)=u(/,?)=0 (2) 

and if uo(x), ui (x) are initial displacement and velocity 

u(jc,0) =Uo6t), ^(x,0) =ui(x), 0<x<l. (3) 
at 

The existence and uniqueness in the scalar case are discussed in [ 12 1 and 1 15 1. 

Arosio and Spagnola |3|, Gough |7| and several other authors referred to in Nayfeh 
and Mook 1 14 1 have studied the problem (la, b) from mathematical and physical points of 
view with a focus on modal analysis and eigenvalues of the system. In particular, Anand 
l2l has investigated (lb) proving that the vibrations are stable even in the planar case. 
The elaborate method of treatment employs stability analysis of the time-dependent Hill's 
equation satisfied by the modal amplitudes. In recent years, Shahruz 1 16 1 has treated the 
problem (la) with a bounded input disturbance and his result shows boundedness of the 
ouput. The approach is systems and control theoretic. 

The mathematical theory of stabilisation of distributed parameter systems is currently 
of interest in view of application to vibration control of various structural elements. The 
most common class of vibration control mechanism is of passive type that absorbs vibra- 
tion energy. A system is called strongly stable if the total energy E(t) of each solution 
of the system converges to zero as time t — > + 00. If the convergence is uniform for t > 
with respect to all initial data in the energy space for which £(0) < °°, the system is called 
uniformly stable. Here our investigation is a direct Lyapunov stability approach to obtain 
uniform exponential energy decay estimate for the system (lb), (2), (3) and establish that 
the system is a passive energy absorber. Such estimate has earlier been obtained by Gorain 
|4| for damped linear wave equation in a bounded domain in R". Gorain and Bose 15 161 
have obtained such estimates for an internally damped beam governed by linear equations 
in torsional and flexural modes of vibration. 



2. Energy of the system 

The total energy at time f, of the system (lb)— (3) is defined by the functional 
1 
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Diferentiating (4) with respect to t and using the governing equation (lb) we obtain 

dx 
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The value of the first integral vanishes in view of the boundary conditions (2). We thus get 



dE 
~d7 



-25 
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du 


Jo 


dt 



dx, f >0. 



(5) 



The negativity of the right hand side of (5) shows that the energy E(t) is decreasing with 
time due to the incorporation of passive viscous damping coefficient 5 > 0. The system 
is thus non-energy conserving. Naturally the question arises as to whether the system is 
uniformly stable or not. An affirmative answer is contained in the following theorem, 
which essentially states that in fact the energy of vibration has uniform exponential decay. 

Theorem. Let u(x,t) be a solution of the system (lb)— (3) with uq € H l [0,l] and Ui G 
L 2 [0, 1], where H l [0, 1] is the Sobolev space of order 1, then the energy E(t) of the system 
defined by (4) satisfies 



£(t)<Me^'£(0), f>0 

for some reals pi > and M > 1. 
From eq. (4) 
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(7) 



From eq. (6) it follows that if E(0) <°° with respect to all initial values uo G H l [0, 1] and 
Ui € L 2 [0, 1] then E{t ) — > as t — > °° and the system is uniformly stable. In order to prove 
(6), let £ > be a fixed constant. Proceeding as in Gorain |4| (see also Komornik II 11 1, 
we define an energy-like Lyapunov functional V according to 



where 



V(t)=E(t) + eG(t), 



G{t) = J^ (u~ + 8\u\ 2 )dx. 



Now since u(0,f) = = u(7,f), the Poincare-like Scheefer's inequality 
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Jo 



|u| 2 dx<^ 
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(8) 



(9) 



(10) 



holds (see |13| and |8|). A simple proof is provided in the Appendix. Hence using 
Schwarz's inequality, we obtain 
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where (10) and the defining equation (4) are used. Also, similarly 



f l , 8l 2 r l 
0<8 |u| 2 dx< / a 

Jo n z a z Jo 



du 



dx 



28l 2 , i 
dx< —.-.Eft), t>0. 



Thus using (11) and (12) in eq. (9) 
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\G(t)\< — [l + —)E(t) = HoE(t) (say), t > 0. 



%a ) 



Also from (9), 



f l du f 
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The inequalities (13) and (14) yield for V (t) (defined by (8)) the estimates 
£/\ 

l-—\E(t)<V(t)<(l + eno)E(t), f>0 



(12) 



(13) 



(14) 



(15) 



where we assume that £ < Tta/l, so that the left hand side of (15) is positive. 

Next, differentiating (9) with respect to f, and using the governing equation (lb), we 
obtain 
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Integration by parts of the first term yields 
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(16) 



If we now differentiate (8) with respect to t and then insert the results of (5) and (16), we 
get 



dV r l 



<9u 



dt 



dx-2eE(t), f>0 



and so if we choose £ < min{ 5, na/l}, 

f<-2eE( t )<--^ 
dt ~ 1 + ejuo 

where (15) has been used. Defining 

2£ 



V(t), t>0 



1 + £jU ' 



(17) 
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the above inequality becomes 
dV 

-j- + nV(f)<0, t>0 
and so multiplying by exp(jiu) and integrating from to t, we obtain 

V(t)<e~^V(0), ?>0. 
Finally, using (15) again, we obtain the estimate (6) where 

M=-i±^L>l. (18) 
1 — el/na 

Hence the theorem. 

Remark. In the above we have also obtained exponential energy decay estimate for the 
solution of the nonlinear damped elastic string for initial conditions uo €//'[(),/], UiG 
L 2 [0,Z]. The decay rate jj, is a function of e as given in eq. (17), where e is the lesser of 
the two quantities 8 and na/l. Since djx/de = 2/(1 + e/io) 2 > 0, fi is maximum when e 
is maximum. In actual cases 8 is the smaller parameter unless the string is very long. So 
from the restriction on £, we can assume e < 8 < na/l. Hence e max = 8 and so 

28 28 

Mmax ~ T+8^ ~ T+TO+M)' (19) 

The right hand side is a function of 81/ na < 1 and is maximum when 81/ Ka = l/y/2. 
Hence 

2 Tta Tta 
^< YT ^= T =0.52 T (20) 

for all possible cases of damping. 

It is of some interest to note that Anand Q] in his modal solution of the problem 
assumed fi = 8. This assumption holds when 8 < na/l. In the contrary case n = na/l 
will suffice exponential decay of the solution. 

The quantity M given by eq. (18) however also increases with £, because dM/de = 
(po + l/na)/{\ - el/na) 2 > 0. Thus for e max = 8, 

M= i±l£±!) = i±M =93 (2D 

\-8l/na y/2-1 ( ] 



when 8l/na= 1/ y2. Thus if we look for steeply falling estimate, the energy gets multi- 
plied many times in the initial stages. 



3. Quasi-steady amplitude estimate 

Scheefer's inequality (10) when applied to the energy functional E(t) given by equation 
(4) yields 
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With the aid of the estimate (6), we obtain 



<M£(0)e 



-Ht 



ft 

Jo 

whose solution yields the estimate for the total amplitude of the string 
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(22) 



We say that the estimate (22) is really close for the quasi-steady stage t > when the 
velocities die down. This is because the kinetic energy part in the derivation of (22) gets 
ignored. Apparantly when t — > +°°, |u| — > 0. 



4. Conclusion 

Here we have established uniform stability of a transversely vibrating string (or a metallic 
wire) fixed at both ends, governed by vector Kirchoff type nonlinear integro-differential 
equation, which takes into account the string's elasticity and passive viscous damping. 
It has been proved directly from the equations of motion that the energy of the system 
decays exponentially with time. Study of nonlinear vibrations assumes significance in 
analysing slender structural elements capable of withstanding finite deformations. This 
paper is motivated by such considerations. 

Appendix 

Here we give a simple proof of Scheefer's inequality (10). Since u(0,f ) = = u(/,r), we 
can have a Fourier representation 

/ \ . / x • nnx 
u(x,t) = 2_ b„(f)sin— . 

n=\ 1 
Assuming term-by-term differentiation holds, 

du % y~i , . nnx 
_ = 7 £ nb „(0cos — . 

Applying Parseval's theorem to the two series, we get 

L 1 n=\ Z 1 n=l 1 J0 

Hence the inequality. 
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